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imnoDucnoN 


An  Electiodieimd-Cliemical  ^TC)  Gun  Modeling  Woikshop  was  held  at  the  Weapons  Technology 
Diiectoiate  (WTD),  Amiy  Reseaidi  Laboratory  (ARLXfonneily  the  Ballistic  Research  Laboratory 
(BRL))  on  12-13  May  1993.  The  wmlcshop  was  spe^cally  focused  on  technical  aspects  of  ETC  gun 
rnoddtaig  and  was  sponsored  the  WTD/ARL.  Thus,  foe  meeting  was  not  Ktuoad  in  sct^  as  an  earlier 
meeting  held  at  BRL  under  the  sponsoflfo4)  of  the  JANNAF  Combustion  Subcommittee  on  9-11  July 
1991. 

The  ot^ecdves  of  foe  workshop  were  to:  (a)  assess  foe  cunent  state  of  ETC  modeling:  (b)  determine  the 
class  of  problems  cuirantly  bei^  addressed  by  various  organizations:  (c)  encourage  exdiange  of  ideas  and 
modeling  techniques  among  participants:  and  (d)  provide  recommendations  for  future  work,  particularly 
in  view  of  decHiting  research  fonds. 

The  following  compUation  contains:  (1)  ttie  meetit^  agenda:  (2)  a  list  of  cunem  known  U.S.  ETC 
modding  activities:  (3)  diagnostics  measurements  desired  by  ETC  modelers  as  summarized  by  the  grotq): 
(4)  presentations  by  organizations  attending  foe  workshop;  and  (S)  a  list  of  participants.  The  workshop 
bcMfitted  greatly  from  tiie  partidpation  of  univetsity,  i^ustry,  and  government  (Army,  DMA.  Navy) 
personnel  and  contractors  all  worked  toward  tte  common  goal  of  understanding  ETC  processes 
through  modeling. 

We  would  like  to  dunk  Ms.  Sharon  Ridurdson  for  her  excellent  administeting  of  the  workshop;  each 
participant  for  dieir  active  roles  during  die  wmkshop;  each  author  and  funding  agency  fwallowirig  work 
to  be  igtrinted  in  ddscompilatioo  (note  that  there  are  undassjfiedAmlimited  and  unclassified^jpvemment 
&  contractors  only  secthms)— and  Army,  IMA  and  Navy  personnel  in  many  difibrem  agencies  for  dieir 
cooperation  and  sifoprut  in  diis  endeavw.  The  modding  commuidty  has  made  significant  advances  in 
understanding  ETC  gun  inaeriwbdlistics  over  the  past  several  years  and  continues  to  impact  gun  cartridge 
design.  It  is  hoped  that  dds  compilation  will  cnivey  die  devetopment  of  statenif-die-ait  models  for  ETC 
guns. 


Ghuia  Wren 
BinOberle 

We^xms  Tectanology  Directorde 
Army  Research  Laboratory 
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Electrodieiinal-Cheniical  (ETQ  Gun  ModeUug  Woikshop  AGENDA 
Weapons  Teclaxdogy  Directed,  U.S.  Anny  Research  Laborstmy 


Wednesday,  May  12. 1993 

8:15  Adminiatnaive  Remarks 

Sharon  Rkhardson 

8:30  Sdenoe  Applications  bnemahonal  CoqxMation,  San  Di^,  CA 

C*C.  lUao,  Fled  Su 

9:30  General  Dynraics  Land  System/GT-Devices 

Nids^^Hnsor 

10:30  Break 

10:45  Sdenoe  ^iplicstions  imentatknal  CorpcMation.  Ft  Washington.  PA 

Sandy  Dadi 

11:45  Lunch 

1:00  Ssndia  National  Laboratories,  Albuquerque,  NM 

Steve  Kempka 

2K)0  FMC  Onpontion,  MhmeipOIis.  MN 

Patridt  Jahke 

3:00  Break 

3:15  U.S.  Amiy  Aimament  Research,  Dev^jpmenc,  and  Engineering  Center 

LeeHairis 

4:15  Pennsylvaiiia  State  Univeisily 

F.  Cheung 

5:00  Adjourn 

7:00  No-host  dinner  at  Tidewater  Grille,  Havre  de  Grace 


nmrsday.  May  13 

8:30  U.S.  Anny  Researdi  Latxnaioiy 

Gloria  Wren 

9M)  CMin  Coiporation 

Hu^  MbElroy 

lOKX)  California  Researdi  and  Tcchndogy 

Phfl  Ifookham 

11:00  Group  discusdon  and  wrap-up 

12M  Adjourn 
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Current  ETC  Modeling  Activities  in  the  U.S.  (May  1993) 
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Diagnostic  Measurements  Desired  by 
ETC  Modelers 


7.  Analysis  of  dosed  bomb  experimental  data  via  7.  Can  be  done  (GTD)  7.  ARL  will  share  dosed 

fluid  dynamic  models  bomb  data  with  GTD 


Diagnostic  Measurements  Desired  by 
ETC  Modelers 
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Intentionally  left  blank. 
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ARL  WORKSHOP  ON  ETC  MODELING 
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Work  performed  under:  DAAA15-90-C-1061 
DAAA1 5-92-0-001 7  and  GDLS  IR&D. 
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Ed  Lilliott,  Hugh  McElroy,  Dennis  Worthington 
and  Gloria  Wren  contributed  either  directly  or 
by  helping  with  codes  or  documentation. 
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COORDINATION  W'TH  OLIN 


ETC  ENHANCEMENTS 
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ETC  OUTPUT  DATA 


CSBNBRAL  DYNAMICS 
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TUTORIAL  IN  PROCESS  AT  ARL 


XNOVAKTC  TRIALS 
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DOCUMENTATION  STATUS 


ci»a>  GENERAL  DYNAMICS 

QT  *  P6Vi06S  land  Systems  Division 
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COMPARISON  WITH  EXPERIMENT 


EARLY  BREECH  IGNITION  CAUSES  AXIAL  WAVE 
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FETC  EXPERIMENTS 


BASIC  BALLISTIC  LAUNCH 


inimg 


PLASMA  CHANNEL 


PHYSICAL  IMPLICATIONS 


CHANNEL  WALL 
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BALLISTIC  LAUNCH 


PROJECTILE  LAUNCH 


FETC  CHEMISTRY  MATCHING 


PRODUCTS  MENTIONED 
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MULTI-DIMENSIONAL  SIMULATION 
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ETC  MODELING  WORKSHOP 
U.S.  Army  Research  Laboratory 
Aberdeen  Proving  Ground.  MD  21005-5066 


CRAFT  FINITE-VOLUME  UPWIND/IMPUCIT  CODE 
Brief  Overview  of  Numerics 
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Lagranglan  Solver/Discrete  Particle  Combustion/Erosion  of 

Particles  from  Surface  —  SP-ETC 

Plasma  Upgrades/Unsteady  •••  Use  of  CRAFT  Numerics 


ETC  MODELING  PHILOSOPHY 
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ETC  GUN  MODELING  REQUIREMENTS 
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CURRENT  PROGRAMS  UTILIZING  CRAFT  CODE 
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MICOM  Ducted  CRAFT/DR  •  Ducted  rocket 

Rocket  •  Carbon  particulate  combustion/dense 

_ ; _ loadings _ 

ARL  LPG  CRAFT/LP  •  Liquid  propellant  gun  simulation 

•  Vortex  sheddInq/LES  research 
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JACOBIAN  OF  THE  FLUX  VECTORS 
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UPWIND  NUMERICAL  FLUX  CALCULATION 
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GEOMETRY 
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Figure  1.  Schematic  of  Coupled  Tuo-Dluenslonal  Model  Problen. 


DENSE-GAS  MODELING  CAPABILITIES 
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:  HIGH  PRESSURE  SHOCK  TUBE  STUDIES 


PLASMA  COUPLING  PROCEDURE 
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INFLOW  MACH  NUMBER  HISTORY  FOR  PLASMA 
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TIME  (IN  MS) 


•  interactioiv  Combustion/ Vaporization  of  Plasma 
and  Woridng  Fluid 


TWO-PHASE  MIXTURE  IS  GOVERNED  BY  AMAGAT’S  LAW 
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TWO-PHASE  EQUATIONS  IN  CRAFT 
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HOW  DO  WE  OBTAIN  THE 
DERIVATIVES  OF  THE  PRESSURE  TERM? 
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USING  THIS  RELATIONSHIP  FOR  IT  IS  STRAIGHTFORWARD  TO 
DERIVE  THE  DERIVATIVES  OF  THE  PRESSURE 


SPEED  OF  SOUND  IN  TWO-PHASE  SYSTEM 
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Speed  of  Sound  in  Liquid 


VALIDATION  STUDIES  FOR  GAS-LIQUID  FORMULATION 
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PRESSURE  PROFILES  FOR  A  TWO-PHASE  SHOCK  TUBE 
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PRESSURE  PROFILES  FOR  A  TWO-PHASE  SHOCK  TUBE 

Gas  =  Low  Pressure  Side  and  Liquid  =  Hig^  Pressure  Side 

MacCormKk  Solulloi.  CRAFT  (RoertVD)  Solution 
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PRESSURE  HISTORY  ON  WALL  OF  TWO-PHASE  SHOCK  TUBE 


PRESSURE  HISTORY  FOR  A  RADIAL  TVi'O-PHASE  SHOCK  TUBE 


FREQUENCY  SPECTRA  FOR  A  RADIAL  TWO-PHASE  SHOCK  TUBE 
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TURBULENCE  CHARACTERISTICS  IN  ULLAGE  TUBE 
17.5  eras  FROM  PROJECTILE  BASE 
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COMPARISON  OF  PRESSURE  HISTORY  WITH 
AND  WITHOUT  SUBGRID  TURBULENCE 
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PRELIMINARY  SOLID  PROPELLANT  BURNING  MODEL 
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COMBUSTION  PRODUCT  FLUX;  p  V  =  p^ 


COUPLING  OF  SOLID  PROPELLANT  BURN  MODEL 

INTO  CRAFT 
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DYNAMIC  MOTION  OF  RECESSING  SOLID  PROPELLANT 
ACCOUNTED  FOR 
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(S/H)  IBA 


TIME  (IN  MS)  TIME  (IN  MS) 


FUTURE  DIRECTION 
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LAGRANGIAN  PARTICLE  SOLVER  OPERATIONAL  IN 
CRAFT 
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ETC  Objectives/Issues 
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Demonstration  calculation:  simulation  of  SNL  experiment  ET038 


Modelling  Focus:  Plasma  Nonuniformities 
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SNL  Capillary  Discharge  Experiments 
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ablating  polyethylene  capillary  liner 


Experimental  Setup 
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Capillary  Discharge  Data 
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Mass  of  ablated  polyethylene:  122  mg 


Energy  Transformation  in  an  ETC  Capiiiary 


Motivation  for  2D  Transient  Model 

•  Why  2D? 
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Model  Description 
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Ideal  Gas  Equation  of  State 


Governing  Equations 
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Ohmic  Heating 

(Electric  Fleld)^*Electrlcal  Conductivity 


Temperature  Evaluation 
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Ablation 


-  Assume  dissociation  temperature  Tq,  density  pp,  iatent  heat  y 

-  Solve  cubic  equation  for  u,  given  radiative  flux 


Ablation  Model  Features 
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Turbulence  in  ETC  Capillaries 
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-  volumetric  heat  source 
*  radiation 


Turbulence  in  ETC  Capillaries,  Cent’ 
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Turbulence  In  ETC  Capillaries,  Cent’ 
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Radiation  dominates  wait*generated  turbutence  for  this  example 

-  radiation  and  turbuience  are  equivaient  at  2  eV,  where  radiai 
variations  (due  to  radiation)  are  negiigibie 


Ibrbulence  in  ETC  Capillaries,  Cont’d 

Turbulent  conductivity  (wail-generated  turbulence) 
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Approach:  use  no  turbulence  model,  examine  solution  fieids  for  waves 


Electric  Field  Model 
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dV/dr 


Plasma  Initiation 
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Uniform  Power  Deposition  Caicuiation 
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>  Current  will  agree  with  data:  power  =  voltage  *  current 
-  Resistance  will  agree  with  data:  power  =  voltage^  /R 


Axisymmetric  Model  of  Capillary  and 
Mixing  Chamber 


Capillary  Radius  0.3  cm 

Overall  Length  1 1 .4  cm 
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Distorted  Geometeiy 


Energy  Deposition  Region 
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100 


Temperature  (4.2  eV) 


200  ^8 


Temperature  (3.3  eV) 
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200  fxs 


Pressure 
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100  ^is 

Polyethylene  Density  (10‘^) 


200  fis 

Polyethylene  Density  (10'^) 


power 


Plasma  Initiation,  cont’ 
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•  Aluminum  vaporizes,  provides  new  electrical  path 

•  Lee-More  Model  for  aluminum  electrical  conductivity  -  progress 


Summary 
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>  provides  high  electricai  conductivity  at  eariy  times 
-  convected  out  of  cartridge  by  late  times 


Doppler  Radar  Measurements 
If  Projectile  Kinematics 


FMC  Naval  Systems  Division 
Minneapolis,  Minnesota 
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Electric  Armamente 
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Electric  Armaments 


Convert  Doppler  frequency  to 
projectile  velocity 


Electric  Armaments 
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Janka/ARL  Prat^l  2/93/9 


Electric  Armaments 
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Electric  Armaments 


Etectrlc  Armamenta 
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Electric  Armaments 
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Inverse  Model  -  Formulation 
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Ruction  of  Propellant  Burned 


JankWARL  Prai/S/12/B3/17 


Electric  Armaments 


Janka/ARL  PraaA/12/B3/16 


Electric  Ammmente 


JankWARL  PraWS/12nV19 


Im-EhmONALLY  LEFT  BLANK. 


140 


Enhanced  Propellant  Bum  Rate 
through  Plasma  Erosion 


L.E.  Harris,  D.  Chul,  J.  Prazelskl,  P.  O'Reilly,  O.S.  Downs 
US  Army  Research,  Development,  and  Engineering  Center 
Plcatinny  Arsenal,  NJ 

•  and 

W.F.  Oberle 

US  Army  Research  Laboratory 
Aberdeen  Proving  Ground,  MD 
and 

J.R.  Grelg 
GT-Devices 
Alexandria,  Va 
and 

H.A.  McElroy  and  J.G.  Buzzett 
Olin  Ordnance 
St.  Petersburg,  FL 


ABSTRACT 

In  both  ETC  gun  firings  and  ETC  closed  bomb  experiments  to  date,  it  has  been 
observed  that  the  bum  rates  of  solid  propellants  are  not  altered  by  the 
presence  of  the  ET  plasma  except  through  the  pressure  created  by  the  addition 
of  the  electric  energy.  However,  in  a  recent  series  of  ET  gun  firings,  we 
have  observed  that  the  bum  rate  of  the  propellant  can  be  dramatically 
increased  by  the  erosive  effect  of  the  plasma.  The  propellant  was  JA-2,  and 
for  these  ETC  gun  firings  the  whole  propellant  charge  was  cast  as  a  "single 
perf"  monolithic  grain.  The  central  perforation  was  0.25  inch  diameter  and 
the  outside  diameter  of  Che  grain  was  1.10  inch.  The  propellant  was  fired  in 
a  30mm  (bore  diameter)  gun  in  the  Army's  ET  Gun  Facility  at  GT>Devices.  The 
plasma  generator  was  a  high  press\ire  capillary  discharge  and  the  emerging 
plasma  was  carefully  injected  into  the  central  perforation  in  the  propellant 
grain.  All  ocher  surfaces  of  the  propellant  grain  were  inhibited  using 
silicone  grease.  The  charge  design  is  described  in  detail,  and  results 
obtained  by  simulating  these  firings  using  IBHVG2  are  dlsctissed.  From  these 
simulations,  it  is  concluded  that  the  bum  rate  achieved  in  these  gun  firings 
was  enhanced  significantly  through  the  erosive  effect  of  the  plasma. 
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Enhanced  Propellant  Burn  Rate  Through 
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GT  Devices 
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CBP  EQUATIONS 
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CBP  EQUATIONS 
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CBP  EQUATIONS 
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Table  III  CBP  Calculotion  of  Interior  Ballistics  Time  Profiles 
in  0  120mm  M256  Tank  Gun 


l(MS) 

Zp(M) 

V(^M/S) 

1 — 

f(t) 

.  M(KG/S) 

0.00 

0.000 

0- 

3,424 

0.004 

0 

0.10 

0.002 

42 

3,244 

0.002 

268 

0.20 

0.008 

85 

3,037 

0.007 

536 

0.30 

0.019 

127  . 

2,933 

0.016 

804 

0.40 

0.034 

169 

2,881 

0.028 

1,072 

aso 

0.053 

211 

2,854 

0.043 

1,341 

1.00 

0.211 

423 

2,812 

0.173 

2,681  , 

1.50 

0;476 

634 

2,804 

0.388 

4,022  ^ 

2.00 

0.846 

846 

2,801 

■0.690 

5,362 

2.40 

1.218 

1,015 

2,800 

0.994 

6,435 

2.41 

1.225 

1,018 

2,800 

1.000 

6,453 
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Figure  4.  Projectile  muzzle  velocity  versus  velocity  for  120*  and  140-mm  tank  guns 
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Figure  6.  120-mm  M291  ETC  tank  gun  cartridge 
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Figure  8.  Muzzle  velocity  In  various  air  defense  guns 


Figure  9.  105-mm  ETC  air  defense  cartridge 
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performance  is  not  achievement,  plasma  influenced  grain 
de-consolidalion  at  the  correct  rale  would  be  required  to  meet 
performance. 
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GLOC  «  -6.5  GLOC  ®  -2.1  .  BASE 

LOCATION 
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SCHEMATIC  DIAGRAM 
30mm  PFN 
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IBHVG2  PREDICTIONS  FOR 


4' 
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JA-2  PRESHOT  IBHVG2  RESULTS 
2.4  MSEC  AT  0.5  MJ/KG  EE 


^  o 

IIJ 

CO 

o  2 


o  o 
o  o 
o  o 


r’’—* 


QUARTER  CHQ - HALF  CHARGE 


Firing  Schedule 
As  Fired 


n 

(A 


n 

>1 

e 

X 


e  e  a  a 

e  e  e  e 

A  AXOA 

i3  AP  A 

9  9  9  9 

9  9  9  9 

H  ^4  H  tc 

£4  £4  tl  04 

*0  *0  “O  *0 

•d  •O  *0  TJ 

• 

0  o  o  o 

O  0  O  O 

a 

ssass 

xxxx 

a. 

o 

in 

00  03 

00  00 

n  n  n  n 

o 

OS  OS 

os  os 

os  os  os  os 

H 

+  + 

+  + 

03  03  00  03  CO  CQ 

+  .+  +  + 

03  03 

00  03  03  03  03  03 

•4  »4  (4  t4  (4  t4 

oQ  n  n  00 

<u 

tJ  a 

.4  iJ 

drjf*dfdrdf*df* 

A  A  tO  A 

O 

1  1 

1  1  1  1  1  1 

1  1  1  1  1  1 

till 

in 

>4 

>  to  C^ 

to  Pt  to  fK  p»  pi| 

to  P*  P*  CM  P^  CM 

to  n  CM  p^ 

9 

M  •  • 

•  •  •  • 

^CM  Pt 

M*  in  PM  m  <n  in 

PM  n  p>  in  CO  in 

^  in  pt  Pt 

9 

•o 

>1 

« 

>lP  g 

A 

Ot<w4  Ot 

A 

pa  '•v 

m  m 

in  in  m  in 

in  in  p*.  p* 

9 

sc  i-3  in  in 

in  in  p*  p*  O'*  r* 

1^  P*  P'’  C" 

p»  c^  CM  n 

* 

C  «  ^  C^  C" 

P'  C'*  PI  P1  •  • 

PM  CM  •  •  •  • 

•  •  •  • 

•a 

04  a  w  .  . 

•  •  •  « o  o 

•  *0000 

o  o  o  o 

9 

u 

•H 

JS 

<M 

PSP 

9  a  cpo-^ 

00 

^  iH  C  a  CM  CM 

CM  PM  iH  H  H  H 

to  to  H  fH  CM  to 

PM  CM  PM  CM 

os 

o  9  a  g  •  • 

•  •  •  • 

0«  0«  1-3  '—CM  PM 

PM  PM  *-l  H  t-l  tH 

b  O  f-l  H  CM  O 

CM  CM  PM  CM 

A 

«  p 

» 

a 

Otj3 

X 

h 

P  o> 

X  c 

«  -H  M"  CM 

M*  PM  PM 

^  CM  ^  PM  CM  ^ 

M>  M* 

C  9 

M4 

j:  9 

9  n 

U  ^  iH  H 

CM  0l<  rM  H  rM  t-l 

H  H  H  H  tH  1-4 

CM  h  CM  Ob 

n 

•0 

A 

o 

u 

9  9 
4J  A  -P 
O  g  « 

x;  9  a 
CO  as  ^ 


m  n 
at  ot 

CO  CO 

o  o 


o  n  PI  n  n  p» 
ot  at  ot  ot  ot  ot 

'«  '«  in  to  to  to 
H  H  rC  H 


p>  n  n  n  n  r> 
ot  ot  ot  ot  ot  ot 

O  O  H  rH  CM  r« 
CM  CM  CM  CM  CM  CM 


CO 


H  CM  ri  <«' in  to  00  ot  o  H  CM  n  ^ 


m  to  r«  CO 

rM  H  rM  H 


P  JS 
M-l  Ot 
9  O 

I  i 

g  g 
o  o 

PI  CM 

I  I 

n  a 

*3  K 


9 

.  c 
P  o 

•H  i-i 

«H  e 
03  03 


s 


+  00 
03  00  h) 
hO  »a.3r 


LB  with  LL9  disconnected  from 


DATE  04/30/93  30MN  AXIAl  CARTRIOCC  OATA  lAK 

06:32  AM  KEViSEO  4/93  fOR  OLIN  URfES  OATA 


:  seS:SKR8!3SSS3C{:C3SK 

m  m  . . . . 

!BS2  ;  o-ooooooooo-o®®®**® 

Mt  O  • 

i  ^ 


«  Ml  «• 

X  ac  % 
u  X  S 


1*1 

^  m 

3S  = 

u* 

ils 


SiSli*  M 

siiic!  Is 

mj  m  m  m  m  m  **x 

s!!j|ss|f 

siiiliisi 

•«  •* 

Jllllljll 

«  •  •  «Q  *  •  •  • 


M  § 

S  S  g  8 
5  8-8 

li  I  jT 

i8|s 

slU 

w  w 

Ml  W  Ml  M 

iSiS 


S3Sss5555 

Sccsiscicctc 


838888888 

ssiaissis 

wou«juSuuu 


sisl 

^  3  ^  K 

8  IC  If  S 
^  •  1^  • 


8  :;8 

9  9  9  9 


O  O  2  <8  o  o 

Z  &  S  Z 

«l  «l  S  M  «• 


S  S  8 
8  ft  8 


0  0  0  19 

o  o  o  o 
%  X  X  X 
o  o  o  o 


<  <  <  < 


<<<<<<<< 


<p»  o  ^  ^  tf%  o 


o 


O  9 

8  : 

> i 3S88S8e888SS8S88SS 


3Si::  =  SSi^8*S8;i;:;38833 


;83S3S;83S8S8838SS 


ggKg8»3Sftft383ft8S§8 


S3S88SSSS8888888S8 


SHOT  Notes  on  retainer  system  Retainer 

NUMB  System 


H  H  H  H  H  H 
H  M  H  H  H  H 
HHHHHHHHHHHHHHHHHH 


«  • 
U  O 

A  e 
o  ft 

•p 

c  A 


■o 

« 

** 

« 


■o 

« 

£ 
o 
•  e 

ft  3 

e  A  • 

WiHH 
P  A 
•  A  D*  e 
*0  JQ  C  « 


A 

e 

e 


P 

A 

X) 

O 

c 


■D 

e 

j; 

u 

c 

3 

A 

1^ 

e 

A 

3 

Ai 


A  O  MH 

h  m-p  e 
A  S  e  op 
«  XI  0<H  •  P 
o  o  a 

A  A  <0  A  A 
^£•0  etH 

0,0  px  0,0 

>  m  -P 


*0  H 

0  A 

Xt  0 

0  0 

•  3 
<H  a 
A 

P  A 

p  ja 

A  3  •  • 

JB  -P  A  A 
O  O 

C  <>  A  A  O  • 

•H  p  H  r-t  •  4J  A 

A  aa'O  c  o 

•O  A  A-H  A 


A 

P 

P 

A 

A 


>*iH 

-I  A 
•P  P 
JS  P 

0»A  _ 

•Hjoxi'oeco3ec 

ft  0  Oft’H  ft  (VH'H 
0  CJSJS 

•HOUAAAAA'O 
•a  C  CX)  AJQX)  A  A, 
A^3333333X: 

j:  o  A  A-p-p-p-p-P  A  A-H  s  c  cx:  ft 

0aftft  3  03  A*H-Ho<oa 

3*J  OiA  wx:  C  A 

a  A  A  A'O'O'O'O'O  r-l<H0AA3jSC 
P  P  A  A  A  A  A  P  a»B  PPAA'H 
A  A  A  APx:x;x:-P  A  CAAH3 
XlXl  C  O  O  OH  C  C  '-H  c  c  a  Dv 

3  3H<H  A  C  C  e  A-HH  e  HH  A  c 
PPAAS333IIA  >0>AAP  D>H 
PP  AAA  PAOBPPACP 
PP  A  A  AHHH  A  APTI-H  a  a  B-H 
ppppr-l  Hpp  pppHpH 

AA  N  0*0>0»N  AP  A  A 

A  AHHNBBBNOAHOHHPOA 
B  BHH  0H<HH  oh  B  AHHH  AHP 

HH<<<asp$ap;20HS(9<<<PSoto 

-Hcin«inu>i>o0oiOHc>in^in«or«>eo 


170 


CO 


^  ^  <s 

«-•  “M- 

o\ 

m 

04 

•  « 
^  VO  Tf 

t> 

CO 

, 

VO 

4-H 

CM 

CM 

cs 

00 

4-^ 

•M* 

VO 

> 

fn 

CM 

cn 

CM 

«n  ov 
^  <s 


00 


o  o 

CO  fO 


w-iocscsr- 

^  (TJ 


(1< 


H 

O 

X 

CO 


«n  r*  c\ 


H  H 


tr> 

CM 


o 

CM 


lO 


- - - i  o' 

o  o  o 

CM 

(r>l)3d 


17 


JA2  MONOLITHS 
HALF  CHARGE 


00  v©  cn  o  00  «n 

00  ro  VO  »n  'O 

^  m  <s  fn 


^  OO  00  «0  Tf 

^  <o  vn  ^  »n  m 


fn  o  o  <N  «n  ^ 
cn  cn  VO  cs  vn  fo 


cu, 


H 

g-Nvo-sns 

CO 


172 


EE(KJ) 


H 

D 

Q  S 
X  H  O 
2  CQ 
<  H 
Pi 

Pt  iJ 

O  £ 
«  H  a 
OU  Pk 


(AH  «OHfHO(f-lH«OHQf-iH(Mr'Hn 
•  O  •••••Ci>***w*0**5C* 
.o2ooooo2ooodo^oooh 
a  K  Ot  K  2 

i4 


•o 

>4 

d) 

3 

c 

*0  *0  ^ 

«0  TJ 

u 

•0  u  o 

■o 

0  0 

3 

o>  c  c  « 

0 

c  c 

Xi 

c  h  h  a 

c 

U  P 

M  3  3  0 

p 

3  3 

3  A  A  O.' 

3 

A  A 

0) 

i3  C  1.  0 

A 

u 

^  % 

0) 

•>  0  *0  *0  04 

•0  TJ 

•H 

n  <o  0  0 

O 

0  0 

a 

0  ax  0 

0  0 

^  0 

U  -P  O  0  D> 

O 

0  0  4-4 

01 

0  rH  0  0  p 

0  0 

0  0  0 

c 

0  0  P  0 

•H  j; 

P  P  A 

o 

ae  0 

0  O  3e 

H(MU(nH<nOtfHOeoHco(no(<MHH 
U  •  O  •  CJ  .  •  .u  >  ‘U  •  •  •  <22 

<a;(M<<n<«;(^coio<4;or'4;«oMo>ooo 
H  K  «  <-(  H  rH  H  H  (O  K  (O  m  o  ZZ 

6h  &<  Eh  H 


•p 

•0  4J  -P  '0 

•0  *> 

P 

•O 

0 

$  a  0  -P 

0 

o 

0 

3 

O  3  3  -P  0 

P  3 

3 

c 

V 

axi  >0  0  3 

0  ■O 

•0 

0 

0  V 

0 

3 

o 

o  o  n  u 

P  n 

H 

J9 

n  H  n  o 

VC 

•o 

4J 

0 

0  « 

0 

c  a 

0  0  0^  •> 

X  0 

0 

C 

P  X 

^XX  C  0 

C  X 

X 

3  C 

C  C  C  3  X 

3  c 

c 

CA  3  P  P 

3  3  3  A  B 

A  3 

3 

0 

x:0043Xjc:o3 

O  X 

j: 

X 

0  O  *0  *0 

OOO  XP4J-P-P  O-P 

0 

o 

e  3 » 

»  o 

0  0  0  a  >  0 

O«oOOoo<-i(M0 

33330io3r- 

PaUlH  OiOiClfHCMO-lir'O'O'O'O^N'OC'l 

'O’lneOriO^COOOt 

)  n  (A  <A  O  O  H  ^ 

to 

<s  r--  (>•  o 

CO  o  in  H 

r-  «o  o  CO  GO  o\  in 

t- 

H  H  lO  lO  CO 

nao«oHciH(v(Aor^in(AOininin'4'r4 

H(nOU)t-ltf)r4«H^COCO(DCOO(SH(0 
'»eo<np*^oo^«^«o(nr*r*n(Mvo<'j(n 
H  H  H  H  H  H 

fHr>i(n^iA(or>co(nOH(Mn^in\or>a 


173 


Residue 


Shot  with  Moderator  Tube 
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-  Some  cracking  in  larger  grains 

More  data  will  be  obtained  in  the  at  ieast  six 
remaining  shots 


INTENTIONALLY  LEFT  BLANK. 
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(Presented  in  the  ETC  Modeling  Workshop  at  the  Army 
Research  Laboratory,  Aberdeen  Proving  Ground,  MD) 
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OBJECTIVES 


•  To  evaluate  the  applicability  of  existing  correlations  for 
describing  the  entrainment  and  subsequent  breakup  of 
the  liquid  propellant  droplets  caused  by  the  discharge  of 
a  plasma  jet  from  a  plasma  generating  cartridge. 


•  To  identify  the  key  elements  of  plasma/liquid  propellants 
interaction  that  are  not  fiilly  understood. 


•  To  summarize  the  findings  of  this  critical  review  and  to 
define  basic  research  needs  for  future  modeling  and 
diagnostic  activities  in  this  subject  area. 


•  To  perform  a  theoretical  analysis  to  provide  guidance  in 
the  development  of  a  suitable  correlation  applicable  to 
ETC  gun  conditions. 
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CHARACTERISTIC  LENGTH  AND 
VELOCITY  SCALES 


Two  Macro-Length  Scales 

dh  —  Hydraulic  Diameter  of  the  Core  Region 
(Associated  with  phenomena  of  the 
Interfacial  Shear  Action  and  Inertia  of  the 
Core  Flow) 

6  —  Liquid  Film  Thickness 

(A^ociated  with  Phenomena  of  Liquid- 
Film  Motion  and  Induced  Surface  Wave 
Form) 

One  Micro-Length  Scale 

/tw  —  Taylor  Wavelength 

(On  the  same  Order  as  Ligament  or  Initial 
Droplet  Size  Before  Secondary  Breakup) 


a 


1/2 


s  (P/  -  Pf) 


Two  Independent  Velocity  Scales 

(  Vf  —  Bulk  Liquid  Film  Velocity 

I  Vr  —  Relative  Velocity  Between  Gas  and  Liquid 
Phases  (  V,  =  Vg  -  ) 

or 


Vf  and  Vg 
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OBSERVATIONS  OF  ENTRAINMENT  PHENOMENA 


•  Similarities  exist  between  the  entrainment  of  liquid 
propellant  by  a  plasma  jet  in  an  ETC  gun  and 
entrainment  of  droplets  in  two-phase  annular  mist  flow. 

•  The  interface  of  the  two  fluids  i"  inherently  unstable. 

•  When  Vr  >  (V,)  critical  or  Vf  >  (Vf)  critical,  instabilities  set  in 
and  grow  in  the  interfacial  region. 

^  formation  of  wavy  interface  and  large 
amplitude  roll  wave  (so-called  Kelvin- 
Helmholtz  instability) 

—  ligaments  and  droplets  tom  off  from  the 
interface  dynamic  interaction 

•  Critical  Conditions  for  the  Onset  of  Entrainment: 
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OBSERVATIONS  OF  ENTRAINMENT  PHENOMENA 

(continued) 


•  Detailed  studies  of  entrainment  mechanism  have  been 
performed  by  Newitt  et  al.  (1954),  Hanratty  and 
Hershman  (1961),  Chung  and  Murgatroyd  (1965),  and 
Ishii  and  Grolmes  (1975). 


•  In  general,  entrainment  may  take  place  in  a  number  of 
different  ways,  depending  on  the  flow  situation. 


•  Hydrodynamic  and  surface  tension  forces  govern  the 
motion  and  deformation  of  the  wave  crests.  Under 
certain  conditions,  these  forces  lead  to  an  extreme 
deformation  of  the  interface  that  results  in  the  breakup 
of  a  portion  of  a  wave  into  liquid  droplets.  The  forces 
acting  on  the  wave  crests  depend  on  the  flow  pattern 
around  them  as  well  as  on  the  shape  of  the  interface. 
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RELEVANT  ARTICLES  ON  ENTRAINMENT 
MECHANISM  STUDIES 


•  "liquid  Entrainment:  The  Mechanism  of  Drop 
Formation  from  Gas  or  Vapor  Bubbles,"  Newitt,  D.  M., 
Dombrowski,  N.  and  Knelman,  F.  H.,  Tran,  Inst,  Chem, 
Engng,,  VoL  32,  p.  244, 1954. 


•  "Fundamentals  of  the  Hydrodynamic  Mechanism  of 
Splitting  in  Dispersion  Process,"  Hinze,  J.  0,,AIChE 
Vol.  1,  p.  289, 1955. 


•  "Initiation  of  Roll  Wave,"  Hanratty,  T.  J.  and  Hershman, 
A,,AIChEJ,,  Vol.  7,  p.  488, 1961. 


•  "Studies  of  the  Mechanism  of  Roll  Wave  Formation  on 
Thin  Liquid  Films,"  Chung,  H.  S.  and  Murgatroyd,  W., 
Symp.  on  Two-Phase  Flow,  Vol.  2,  Paper  A2,  Exeter,  UK, 
1965 


•  "Prediction  of  Onset  of  Entrainment  for  Liquid  Metals," 
Ishii,  M.  and  Grolmes,  M.,  Trans,  Am,  Nucl  Soc,,  Vol. 
21,  p.  325,  1975 
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Type  I 
Roll  Wave 


Type  2 

Wave  Undercut 


'//////////////////////, 


Types 

Bubble  Burst 


'//////////////////////, 


Type  4 

Liquid  Impingement 


Types 

Liquid  Bulge  Disintegration 
(Counter  Current) 


A  Schematic  Diagram  Showing  the  Rve  Types  of  Entrainment  Mechanisms 
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FIVE  BASIC  TYPES  OF  ENTRAINMENT 
MECHANISMS  FOR  TWO-PHASE  ANNULAR  MIST 
FLOW  [ISHH  AND  GROLMES  (1975)] 

ROLL  WAVE: 

The  drag  force  acting  on  the  wave  tops  deforms  the 
interface  against  the  retaining  force  of  the  liquid 
surface  tension.  The  tops  of  the  large  amplitude 
roll  waves  are  sheared  off  from  the  wave  crest  by 
the  gas  flow  and  then  broken  into  small  droplets, 
(most  relevant  to  ETC  gun  conditions) 

WAVE  UNDERCUT: 

Entrainment  is  caused  by  the  undercutting  of  the 
liquid  film  by  the  gas  flow.  This  mechanism  is 
similar  to  droplet  disintegration  by  the  gas  stream, 
(could  also  happen  in  ETC  gun  condition) 

BUBBLE  BURST: 

Associated  with  the  bursting  of  gas  bubbles.  The 
droplets  may  be  generated  by  the  bubble  rising  to 
the  surface  of  a  liquid. 

LIQUID  IMPINGEMENT: 

This  is  caused  by  the  impingement  of  relatively  large 
liquid  droplets  to  the  mm  interface  for  production 
of  small  droplets. 

LIQUID  BULGE  DISINTEGRATION: 

When  a  counter-current  flow  reaches  the  flooding 
condition,  large  amplitude  waves  can  emerge  from 
the  liquid  film  and  coalesce  to  from  a  bridge. 
Bulged  regions  of  the  bridge  can  then  disintegrate 
into  small  droplets  due  to  the  gas  dynamics,  (not 
relevant  to  ETC  gun  conditions) 
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ONSET  OF  KELVIN-HELMHOLTZ  INSTABILITY 


When  both  the  gas  flow  and  the  induced  liquid  film 
motion  are  sufficiently  high,  the  interfacial  waves 
transform  into  large  amplitude  roll  waves. 

Beyond  this  critical  point,  the  interfacial  shear  forces 
become  greater  than  the  surface  tension  forces,  and  the 
onset  of  entrainment  occurs. 


LIQUID  FILM 


ROLL  WAVE 


ONSET  OF  ENTRAINMENT 


•  The  criterion  developed  by  Ishii  and  Grolmes  (1975)  was 
based  upon  a  force  balance  at  the  crest  of  roll  waves. 
Droplet  entrainment  starts  when 

Fi  >  F, 

interfacial  surface  tension 

.  shear  force  L  retaining  force  . 

•  The  dimensionless  critical  superficial  gas  velocij^  at  the  onset 
of  entrainment  was  found  to  depend  on  the  film  Reynolds 
number  and  the  viscosity  number  of  the  fluids,  i.  e. 


(Jg)  critical  ^  (F.ef,  N^) 


where  Pf  and  \Xf  are  the  density  and  viscosity  of  the  liquid,  du 
the  hydraulic  diameter  of  the  cone,  and  if  the  superficial 


•  Jg  is  defined  as 


ONSET  OF  ENTRAINMENT  (continued) 

•  For  a  horizontal  two-phase  annular  flow  with 

160  <  Ref  <  1635,  the  criterion  of  Ishii  and  Grolmes 
(1975)  gives 


4  i  11.78  Rej'!^ 

Jg  i  1.35 

fo'K.  >  ^ 

•  For  values  of  Ref  that  are  above  1635,  a  so-called 
completely  rough  turbulent  regime  occurs  in  which  the 
critical  superficial  gas  velocity  becomes  independent  of 
the  film  Reynolds  number.  According  to  Ishii  and 
Grolmes  (1975),  the  inception  criterion  for  this  flow 
regime  is  given  by 


00 

Al 

^  0.1146 

•  The  inception  criterion  compares  favorably  with  a  large 
number  of  experimental  data  for  various  ^es  of  fluids 
covering  a  wide  range  of  the  film  Reynolds  numbers  and 
the  viscosity  number. 

•  The  inception  criterion  of  Ishii  and  Grolmes  should  be 
applicable  to  the  case  of  plasma/liquid  propellants 
interaction  provided  that  fluid  properties  are  known. 
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Comparison  of  the  Inception  Criterion  of  Ishii  and  Grolmes  with  Data 
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METHODS  FOR  MEASURING  THE  FRACTION 
OF  LIQUID  MASS  FLUX  ENTRAINED 
IN  THE  GAS  CORE 


•  Local  Probe  Measurements 

—  To  determine  the  axial  liquid  mass  flux  at  the 
location  of  the  sampling  probe. 

—  Used  with  limited  success  by  Wicks  and  Dukler 
(1960),  Magiros  and  Dukler  (1961),  Wailis  (1962), 
Steen  and  Wallis  (1964),  Cousins  et  al.  (1965),  Gill 
and  Hewitt  (1968),  and  Yablonik  and  Khaimov 
(1972). 

—  Measurement  along  the  centerline  with  the 
assumption  that  the  mass  flux  is  radially  uniform. 


•  Liquid  Film  Row  Measurements 

—  Row  rate  determined  by  removing  liquid  film 
completely  from  the  test  section. 

—  This  method,  which  eliminates  those  uncertainties 
associated  with  the  local  probe  measurement,  is 
probably  a  more  accurate  measurement  technique. 

—  Successfully  employed  by  a  number  of  investigators 
including  Paleev  and  co-workers  (1962,  1966), 
Cousins  and  Hewitt  (1968),  Petrovichev  et  al. 
(1971),  and  Ishii  and  Mishima  (1981). 
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ENTRAINMENT  CORRELATIONS 


•  Purely  Empirical  Approach 

—  Wicks  and  Dukler  (1960) 

—  Minh  and  Huyghe  (1965) 

—  Paleev  and  Filipovich  (1966) 

-  WalUs  (1968) 

•  Semi’Empirical  Approach  Based  Upon 
Mechanistic  Models 

-  Hutchinson  and  Whalley  (1973) 

—  Dallman,  Laurinat,  and  Hanratty  (1984) 

-  Ishii  and  Mishima  (1989) 
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ENTRAINMENT  CORRELATION  OF  WICKS 
AND  DUKLER  (1960) 


•  Developed  a  correlation  that  relates  the  entrainment 
parameter,  R,  to  the  Martinelli  parameter  X,  i.  e. 

R  =  R(X) 


•  Entrainment  Parameter: 


idPldz\ 


where  is  the  liquid  drop  mass  flow  rate,  and  Weg  the 
critical  Weber  number  defined  in  terms  of  the  film  thiclmess 
6  by 


where  Vg  is  the  velocity  of  the  gas  phase. 


•  Martinelli  Parameter 


(dPldz\ 

j3m% , 


in 


where  quantities  in  the  numerator  and  denominator  represent 
the  sinme-phase  pressure  drops  which  would  exist  if  each 
•  phase  flowed  alone. 

•  Two  Major  Drawbacks 

—  The  dependence  of  the  entrainment  on  various 
controlling  factors  is  hidden  by  the  use  of  the  Martinelli 
parameter. 

—  The  correlation  is  dimensional  which  severely  limits  the 
range  of  applicability. 
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ENTRAINMENT  CORRELATION  OF  MINH 
AND  HUYGHE  (1965) 


Empirical  correlation  relates  the  entrained  fraction,  £,  to  the 
gas  core  momentum,  i.  e. 

£=^(pw7) 


where  Pg  is  the  homogeneous  density  of  the  gas  core  given  by 


Pr  =  Pr 


*  1 
P*^*  J 


where  jd  is  the  superficial  velocity  of  the  entrained  liquid 
droplets. 


•  Entrained  Fraction: 


where  is  the  mass  flow  rate  of  the  droplets  and  Wf  the 
total  liquid  mass  flow  rate  including  both  the  liquid  film  and 
the  droplets. 

•  Two  Major  Drawbacks 

~  The  use  of  the  dimensional  parameters  results  in 
uncertainties  regarding  the  dependence  of  the 
correlation  on  the  fluid  properties. 

—  The  range  of  applicability  is  virtually  unknown  as  the 
effects  of  many  important  factors  such  as  the  liquid 
Reynolds  number  are  not  included. 
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ENTRAINMENT  CORRELATION  OF  PALEEV 
AND  FIUPOVICH  (1966) 


•  The  shortcoming  in  the  Minh  and  Huyghe  correlation 
(1965)  was  eliminated  by  introducing  a  dimensionless  gas 
flux*  in  correlating  the  entrained  fraction  in  the  following 
form. 


=  0.985 


-  0.44  log 


p 

[Hr  41 

2 

X  10" 

.  p/ 

1  O  j 

■ 

where  W®  is  the  mass  flow  rate  of  the  film  and  Wf  the 
total  liquid  flow  rate  including  the  film  and  droplets. 


•  Entrained  Fraction: 


E 


•  This  correlation  showed  fairly  good  agreement  with  a 
limited  number  of  data. 

•  Major  Drawback 

—  The  range  of  applicability  is  not  well  defined  as  the 
important  effects  of  the  hydraulic  diameter  and  the 
Reynolds  number  were  not  included. 


The  dimensionless  gas  flux  is  very  similar  to  Jg  used  in 
the  onset  of  entrainment. 
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ENTRAINMENT  CORRELATION  OF  WALLIS  (1968) 


•  In  an  effort  to  improve  the  correlation  of  Paleev  and 
Filipovich  (1966),  Wallis  introduced  a  modified 
dimensionless  gas  flux  (tz)  by  replacing  the  liquid 
viscosity  by  the  gas  viscosity  in  the  Paleev  and  Filipovich 
correlation,  i.  e. 


E  =  E  (ti) 

where 


n  = 


f  V  u  'I 

8 

'p*  I 

1  o  J 

[p/  J 

•  Major  Drawbacks 

—  The  range  of  applicability  is  not  well  defined. 

--  The  dimensionless  gas  flux  is  not  general  enough  to 
account  for  the  effect  of  fluid  properties. 
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ENTRAINMENT  CORRELATION  OF  HUTCHINSON 

AND  WHALLEY  (1973) 

•  Proposed  a  mechanistic  model  to  correlate  the  rates  of 
entrainment  and  deposition  of  droplets  to  the  liquid  film. 

•  The  deposition  rate  is  taken  to  be  a  linear  function  of  the  droplet 
concentration  in  the  gas  core,  i.  e., 

b  ^kC 


where  k  is  the  mass  transfer  coefficient. 

•  The  entrainment  rate  is  considered  proportional  to  the  equilibrium 

concentration,  i.  e. 

E^kC, 


The  equilibrium  concentration  is  correlated  in  a  functional  form 
given  oy 

C  =  C 


where  t;  is  the  interfacial  shear  of  the  two-phase  flow  and  5  the 
liquid  film  thickness. 

•  Major  Drawbacks 

—  It  is  difficult  to  accurately  determine  the  equilibrium 

concentration.  The  data  for  this  quantity  showed 
considerable  scattering,  sometimes  up  to  more  than  an  order 
of  magnitude. 

—  In  order  to  estimate  the  entrainment,  two  different 

correlations,  i.  e.,  one  for  the  entrainment  rate  and  the  other 
for  the  deposition  rate,  need  to  be  employed.  Any  error  in 
determining  the  deposition  would  add  to  the  error  associated 
with  the  entrainment  rate. 

—  Conceptually,  the  entrainment  rate  is  not  dictated  by  the 
equilibrium  concentration. 
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ENTRAINMENT  CORRELATION  OF  DALLMAN, 
LAURINAT  AND  HANRATTY  (1984) 

•  Conducted  an  experimental  study  of  horizontal  annular  two- 
phase  flow  using  air  and  water. 

•  The  amount  of  liquid  entrained  in  the  gas  flow  is  assumed  to 
be  governed  by  a  dynamic  balance  oetween  the  rate  of 
atomization  of  the  liquid  film  and  the  rate  of  deposition  of 
droplets  from  the  gas  core. 

•  The  rate  of  deposition  of  liquid  droplets  was  considered  to 
vaiy  linearly  with  the  concentration  of  the  droplets  whereas 
the  rate  of  atomization  of  the  liquid  film  was  assumed  to  vary 
linearly  with  the  liquid  film  flow  rate. 

•  Entrainment  Correlation  Developed: 


E  ^  -2&) 


where  A  is  a  dimensional  constant  having  the  value  of 
3.6  X  10-^  in  SI  unit  and  the  maximum  entrainment 
fraction.  The  latter  quantity  is  ^en  by 

£  =  1  - 


where  (Wff)c  is  the  so-called  critical  film  flow  rate. 

•  Critical  Film  Rowrate 

—  The  critical  film  flow  rate  has  been  determined  by 
Dallman  et  al.  (1984)  as  a  function  of  the  ratio  between 
the  droplet  mass  flow  rate  and  the  mass  flow  rate  of  the 
gas. 


206 


Variation  of  the  Critical  Film  Flow  Rate  with  the  Mass-Flow-Rate 
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ENTRAINMENT  CORRELATION  OF  DALLMAN, 
LAURINAT  AND  HANRATTY  (1984) 
(continued) 


•  Major  Drawbacks: 

—  The  correlation  is  given  in  dimensional  form;  thus,  the 
applicability  is  limited. 

—  The  effect  of  fluid  properties  are  not  included. 

—  While  the  concept  of  critical  film  flow  rate  is  correct,  the 
functional  dependence  of  (Ws)c  is  generally  unknown. 

—  The  correlation  is  cumbersome  as  and  8  are 
unknown  quantities. 
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ENTRAINMENT  CORRELATION  OF  ISHII  AND  MISHIMA 

(1989) 


•  Based  upon  the  mechanism  of  roll  wave  entrainment:  The 
shearing  off  of  roll  wave  crests  was  considered  to  be  the 
dominant  mechanism  of  liquid  entrainment. 

•  Droplet  entrainment  would  result  whenever  the  retaining 
force  of  surface  tension  is  exceeded  by  the  interfacial  shear 
force. 

•  The  same  forces  controlling  the  onset  of  entrainment  were 
assumed  to  control  the  entrainment  itself.  However,  a 
modification  of  the  inertia  of  the  gas  core  flow  was  made  to 
account  for  the  effect  of  droplet  inertia,  since  there  are  many 
droplets  flowing  in  the  gas  core. 

•  Under  an  equilibrium  condition  in  which  the  deposition  rate 
exactly  balances  the  entrainment  rate,  the  entramed  fraction 
was  correlated  in  the  following  form; 

E  =  tanh  (7.25  x  lO*’  We'-^  Re,-^) 


where  We  is  an  effective  gas-phase  Weber  number  defined  by 


We 


•  The  Ishii  and  Mishima  correlation,  which  compares  favorably 
with  a  large  number  of  data  ,  is  probably  the  best  correlation 
available  today. 

•  Unfortunately,  the  entrained  fraction  represents  the 
equilibrium  value  that  includes  the  contribution  due  to 
droplet  deposition.  Thus,  this  correlation  tends  to 
unaerestimate  the  actual  entrainment  rate  in  the  ETC  gun 
environment. 
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ENTRAINED  FRACTION 
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Comparison  of  the  Entrainmoit  Correlation  of  Ishii  and  Mishinia  with  Data 
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SPECIAL  FEATURES  OF  DROPLET  ENTRAINMENT 
IN  AN  ETC  GUN  ENVIRONMENT 

•  Droplet  life  times  are  relatively  short  due  to 

—  High  rates  of  heat  transfer 

—  Evaporation,  thermal  decomposition,  and 
combustion 

—  Secondary  breakup  to  smaller  droplets  by  highly 
turbulent  flow  conditions 

•  The  deposition  rate  is  negligible  since  the  droplets  have 
very  litde  chance  to  return  to  liquid  film  before  they  are 
totally  consumed.  The  absence  of  the  deposition  process 
in  an  ETC  gun  represents  a  unique  feature  which  is 
different  from  all  aimular  two-phase  flows  studied  by 
previous  researchers. 

•  For  ETC  guns  utilizing  piccolo  igniters,  the  phenomena 
of  mixing  of  plasma  and  liquid  propellants  involve  two 
sequential  stages. 

1 )  Development  of  numerous  mini-Taylor  cavities 
and  subsequent  expansion  of  cavity  sizes. 

2)  Merging  of  mini-Taylor  cavities  onto  a  single 
large  cavity. 
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RESEARCH  NEEDS  ON  DROPLET  ENTRAINMENT 


•  The  single  most  important  element  in  accurate 
prediction  of  the  performance  of  ETC  guns  is  the 
droplet  entrainment  rate,  since  it  controls  the  rate  of 
conversion  of  chemical  energy  into  propulsive  energy. 

•  A  suitable  entrainment  correlation,  excluding  the  effect 
of  droplet  deposition,  should  be  developed  for  ETC  gun 
application.  This  can  be  accomplished  by  performing  an 
experimental  study  to  measure  the  rate  of  entrainment 
for  the  case  without  the  effect  of  deposition. 

•  In  order  to  develop  a  realistic  model  for  simulating  the 
mixing  phenomena  of  plasma  jets  and  liquid  propellants 
in  ETC  guns  using  piccolo  igniters,  it  is  useful  to  conduct 
detailed  observations  of  the  evolution  of  mini-Taylor 
cavities. 
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THEORETICAL  ANALYSIS  OF  DROPLET 
ENTRAINMENT  IN  AN  ETC  GUN  ENVIRONMENT 


•  The  starting  point  of  the  analysis  is  the  selection  of  an 
appropriate  inception  criterion.  Since  the  roll  wave 
mechanism  is  the  prevailing  mechanism  of  entrainment 
during  plasma/liquid  propellant  interaction  in  an  ETC 
gun,  it  is  deemed  appropriate  to  adopt  the  inception 
criterion  of  Ishii  and  Grolmes  (1975)  in  the  present 
analysis. 

•  Owing  to  the  relatively  high  velocity  of  the  plasma  jet, 
the  plasma/liquid  propellant  system  in  an  ETC  gun  is 
likely  to  be  significantly  above  the  inception  criterion. 
However,  there  is  a  potential  for  the  liquid  film  Reynolds 
number  to  be  reduced  to  a  critical  value  as  the  liquid 
film  thickness  decreases.  Below  this  critical  RCf,  no 
further  entrainment  is  possible. 
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THEORETICAL  ANALYSIS  OF  DROPLET 
ENTRAINMENT  IN  AN  ETC  GUN  ENVIRONMENT 

(continued) 


•  An  upper  theoretical  limit  for  the  entrained  fraction  may 
thus  be  determined  by  assuming  that  the  liquid  film  flow 
rate  is  given  by  the  critical  condition  for  the  onset  of 
entrainment: 


11.78  ATo  ® 


\l/2 


where  (Ref)c  is  now  the  critical  Reynolds  number  on  the 
critical  film  flow  rate  introduced  initially  by  Dallman, 
Laurinat  and  Hanratty  (1984). 

•  The  critical  Reynolds  number  is  defined  by 

(Ee,X  = 


where  (jff)c  is  the  critical  superficial  veloci^  of  the  liquid 
film.  The  use  of  this  definition  in  the  above  equation  is 
equivalent  to  ignoring  the  efrect  of  droplets  in 
determining  the  critical  condition. 


•  In  terms  of  the  liquid  film  Reynolds  number  (Ref),  the 
critical  superficial  velocity  of  the  liquid  film  is  given  by 


(Jffi  =  ^ 


11.78  * 


I 


s 
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THEORETICAL  ANALYSIS  OF  DROPLET 
ENTRAINMENT  IN  AN  ETC  GUN  ENVIRONMENT 

(continued) 

•  Following  the  conventional  definition,  the  entrained 
firaction,  £,  can  be  expressed  in  terms  of  the  superficial 
velocities  of  the  total  liquid  (including  the  liquid  film  and 

the  entrained  droplets)  and  the  liquid  film  alone,  i.  e., 

•  •  •  • 

E  =  ^  =  =  1  - 

Jf  Jf  Jf 

•  The  upper  limit  of  the  entrained  fraction  is  given  by 


After  substituting  the  expression  of  (jff)c  into  the  above 
equation,  Eup  can  be  expressed  as 


11.78 


3 


•  Theoretically,  the  actual  entrainment  during  the  mixing 
of  plasma  and  liquid  propellants  in  an  ETC  gun  should 
be  bounded  by  the  above  upper  limit  and  the  lower  limit 
given  by  Ishii  and  Mishima  (1989),  i.  e. 

(^)jM  <  ^  <  ^up 
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THEORETICAL  ANALYSIS  OF  DROPLET 

ENTRAINMENT  IN  AN  ETC  GUN  ENVIRONMENT 

(continued) 

From  the  above  equations,  a  functional  form  can  be  postulated  for 
the  actual  entrained  fraction  as 


E  =  1 


-  B  Ref'  [ 


\3/2 


/ 


where  B  is  a  coefficient  that  needs  to  be  determined 
experimentally. 

In  view  of  the  fact  that  the  critical  superficial  velocity  of  the  liquid 
film  represents  the  lower  limit  of  the  liquid  film  flow,  the  value  of 
B  should  be  bounded  by 

B  ^  1.635  X  l(f 

In  terms  of  the  dimensionless  characteristic  gas  velocity  as 

Hr 

«  Pf!  P*  > 


a  more  general  expression  for  the  entrained  fraction  can  be 
postulated.  Tliis  is 

E  =  1  -B  N^Re/Jl 


Based  upon  the  preceding  physical  arguments,  the  values  of  the 
exponents  should  be  approximately  equal  to:  a  =  2.4,  b  =  -1,  and 
c  =  -3.  Experimental  data  to  be  obtained  in  the  proposed 
experiments  will  be  used  to  obtain  a  final  expression  for  the 
entrainment  correlation. 


ENTRAINED  DROPLET  SIZE  DISTRIBUTION 


•  In  an  ETC  gun  environment,  the  chemical  reactions  that 

are  directly  responsible  for  the  transfer  of  chemical 
energy  to  propulsive  energy  are  controlled  not  only  by 
the  entrainment  rate  but  also  by  the  entrained  droplet 
size  distribution. 


•  In  an  ETC  gun  chamber,  the  majority  of  droplets  are 
generated  by  entrainment  and  not  by  secondary  breakup 
mechanisms.  Examples  of  secondary  breakup  are 
turbulence  eddy-induced  breakup  and  particle/particle 
collision  generated  breakup. 

•  The  size  distribution  should  be  governed  by  the  process 
of  the  droplet  entrainment  based  on  shearing  on  of  the 
roll  waves. 


Studies  of  droplet  sizes  have  been  performed  by  Hinze 
(1955),  Hass  (1964),  Wicks  and  Dukler  (1966),  Cousins 
and  Hewitt  (1968),  Tatterson  et  al.  (1977),  and  Kataoka, 
Islui,  and  Mishima  (1983). 

A  correlation  of  droplet  size  distribution,  which  could  be 
applied  to  the  interior  ballistic  processes  of  an  ETC  gun, 
was  obtained  by  Kataoka,  Ishii,  and  Mishima  (1983). 
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DROPLET  SIZE  CORRELATION  BY  TATTERSON 
DALLMAN  AND  HANRATTY  (1977) 


•  An  experimental  study  was  performed  to  measure  the 
entrained  droplet  size  m  an  annular  two-phase  flow  using 
a  log-normal  distribution  function  to  simulate  the  droplet 
size  distribution: 


Weid^J  =  0.106 


xl/2 


where  Pg  is  the  gas  viscosity,  Re^  the  gas  Reynolds 
number,  and  'Wt(ayn)  the  charactemtic  V^ber  number 
based  on  the  volume  median  diameter. 

•  The  gas  Reynolds  number  and  the  characteristic  Weber 
number  are  defined  respective^  by 

Q  j  dj.  Q  d 

Re  =  Ilh-JL  and  We  (d,J  = 

«  ^  '  o 


where  dvm  is  the  volume  median  diameter.  Physically,  dvm 
is  defined  such  that  droplets  having  diameters  greater 
than  the  volume  median  diameter  would  occupy  exactly 
half  of  the  total  droplet  volume  in  the  gas  core. 

•  Major  Drawback: 

Their  assumption  of  a  potential  flow  is  not  strictly 
valid  for  two-phase  annular-mist  flow  with  a  highfy 
turbulent  gas  core.  As  a  result,  the  correlation  of 
Tatterson  et  al.  does  not  compare  favorabty  with 
experimental  data  over  a  wide  range  of  gas 
Reynolds  number. 
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DROPLET  SIZE  CORRELATION  BY  KATAOKA 
ISHn  AND  MISHIMA  (1983) 


•  A  mechanistic  model  was  developed  that  was  based  on 
the  mechanism  of  shearing  off  of  the  roll-wave  crest  for 
generation  of  liquid  droplets  by  a  viscous  shear  flow, 
^e  potential  flow  assumption  employed  in  the  standard 
Weber  number  criterion  was  eliminated  in  developing 
the  model  development. 

•  By  incorporating  the  mechanistic  model  with  test  data, 
Ifataoka,  Ishii,  and  Mishima  (1983)  successfully  derived 
a  correlation  for  deducing  the  volume  median  diameter 
from  the  Weber  number 


We{d^J  =  0.028 


1 

-1/3 

•V'l 

p/  J 

J 

2/3 


•  In  terms  of  the  characteristic  Weber  number,  the  volume 

median  diameter  is  given  by 

_  g  We{d,J 

The  above  expression  correlates  the  e^erimental  data 
of  volume  median  diameter  within  ±  40%  errors. 
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Correlation  of  the  Volume  Median  Diameter  by  Kataoka,  Ishti  and  Mishima 
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DROPLET  SIZE  CORRELATION  BY  KATAOKA 
ISHn  AND  MISmMA  (1983) 

(continued) 

The  mean  values  of  the  data  on  drop  size  distribution 
were  fitted  by  Kataoka,  Ishii,  Mishima  (1983)  to  the 
upper  limit  log-normal  distribution: 

^  =  -_L  exp[-(5y)^] 


where  A  is  the  volume  fraction  oversize  defined  as  the 
volume  fraction  of  droplets  having  a  diameter  larger  than 
a  given  diameter,  d. 

In  the  above  equation,  $  is  a  distribution  parameter  and 
y  is  defined  in  terms  of  the  maximum  diameter  and 
volume  mean  diameter  dy^  as 


y 


k  d 

(d:-d) 


and 


Based  upon  the  experimental  results,  the  values  of  £  and 
k  were  found  to  be  5  =  0.884  and  k  =  2.13  whereas 
using  the  upper  limit  log-normal  distribution,  the 
maximum  diameter  was  related  to  the  volume  mean 
diameter  by 


d^  =5  3.13  dym 
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DROPLET  SIZE  CORRELATION  BY  KATAOKA 
ISmi  AND  MISHIMA  (1983) 
(continued) 


•  The  correlation  of  drop  size  distribution  developed  by 
Kataoka,  Ishii,  and  Mishima  (1983)  is  probably  the  most 
suitable  one  among  the  existing  correlations  for  use  in 
simulating  the  phenomena  of  mixing  of  plasma  and 
liquid  propellants  in  the  ETC  gun  environment. 

•  It  is  recommended  that  the  correlation  be  used  in  future 
modelmg.  The  reasons  are: 

—  The  correlation  is  based  upon  a  sound  physical 
model,  taking  full  account  of  key  elements  involved 
in  the  droplet  formation  process 

—  The  use  of  the  upper  limit  log-normal  distribution 
function  is  evidently  appropnate  for  annular-mist 
flow 

—  The  correlation  compares  veiy  well  with  available 
experimental  data  on  drop  size  distribution  covering 
a  wide  range  of  conditions 

—  The  correlation  is  easy  to  use  since  the  primaiy 
quantity  that  needs  to  be  determined  is  the  volume 
median  diameter. 
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Electrothermal-Chemical  Gun 
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Interior  Bailietic  Component  Processes 
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Dimensionality  for  IB  Simulations 
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ETC  IB  Models 


Model 

Description 

Status 

P2SM 

00,  power  syatams  for  electric  guns, 
includes  mission  scenarios 

Compiele 

Piasma 

10  steady  state  capillary 

10  transient  capillary 
•  treats  metals  (In  Progress) 

10  piccolo  tube 

20  capillary 

Complete 

Complete 

In  Progress 

In  Progrm  (SNL) 

BLAKE-ETC 

Thennochendcal  code,  treats 
electrical  energy  as  a  constituent 

Complete 

IBHVQ2-ETC 

OO,  EE  source  term 

Input  power 

-  coupled  to  plasma  (In  Progress) 

Complete 

SPETC 

OO,  SP,  Oetermlnes  power  pulse  for 
given  performance 
-  of^miaadon  (In  Progress) 

Complete 

LUMPET 

00,  LP,  Arrefienius  EOS 

EE  source  term 

Complete 

Conpraaa 

00,  Constant  pressure  solution, 

EE  source  term 

Complete 

XKTC-ETC 

10,  Solid  propellant  code, 

EE  source  term 

Complete 

CRAFT 

10/20/30  research  code,  fundamental 
physics,  state-of-the^rt  numerics 
-  gas/gas  (Complete) 

•  gaa/liquid  (Complete) 

•  gaa/llquid/solid  (In  Progress). 

In  Progress 
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Power  Model 
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Status  of  Plasma  Activities 
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Initiated  calculations  for  breech-fed 
injector  accounting  for  vaporization  of 
A1  fuse 


Experimental  Fixture 


THIN  PLASTIC  FILM 


Ignition  Characteristics 
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Closed  Chamber  Burn  R 


Experiment  0CB14 
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-  Ramped  pulse  over  1.9 
•  41.167  Uj 


M43  Burn  Rate  with  Plasma 


Shots  160  and  161  from  Soreq 
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-  unarge  weignt:  i  ./3  Kg 

-  Propellant:  M30,  7perf 

-  Projectile  mass:  1 .35  Kg 


Power  vs  Time 


(•ui)  *011 X 


#1-  Propellant  burning  rate  (cm/s):  r=.183796P’ 


Simulation  of  OCB1 4 
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( inittol  charge:  6.5.3  g  M43  propeHant) 


Radial  Structure  of  Temperature 
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Intentionally  left  blank. 
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ORDNANCE  imooamBircoitrMmoim 
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ORDNANCE  rmooeumBncimMmomeomKiumoHcom>emioHtBunmmwi»M.9 


WHY  MODEL  -  AT  OLIN 
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comjmmoimctmpoRAnoMcomKnnoiiaBisnmiurtHALM 


WHAT  WE  MODEL  ROUTINELY  AT 
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CONTMW  OUN  OORmMTNlN  OOMIPCnnON  SENSITIVE  aUTEflMit 


LUMPED  PARAMETER 
0-D  IB  Code  •  IBHVG2.504G6 
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OHDNANCE  nBtoocummcomm$ouMC(mH)RAytoMC(m>tmioM$eitsnwe»ukTemM.» 


5  DATA  VS  CODE  RESULTS 
DC  26-PCI 


mcoNTAMSMKcommTKMcattnTinoH-seNsiTiveMAn/m. 


30mm  OLIN  EXP.  DATA  vs.  CODE  RESULTS 


imium  3AUISN3S-NOUU3dHOO  HUlWOdUOO  MW  SmiMOO  IN: 


1-D  XNOVAKTC  (modified) 
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ORDNANCE  rmoocmmiTcomMmoLm 


o 

■sCOMM  This  ’-'jn  converted  to  XNOVdk  TC  3/1^' ^3  K  RL 
N0889.1N  Commented  out  et  stuff. 

10mm  subscale  ARTILLERV  REVERSE  ENGINEER  fMG  FOR  ER' -u.* 
iJSING  ADJUSTED  -^TM  LAvER  TRANSIT  lOM 
ADJUSTING  PATIO  OF  LARGE  TO  SMALL  SPAIN  CICE3  TO  DEMONSTRATE 
PROCESS  CONTROL 

t  IT(F3 

RUN  *  30mm  ARTILLERy  RESULTS 


SXCONTROL 


*M2  CONTROL  DATA 


NPR I NT 

3 

1 

NGRAPw 

=  n 

NOISK 

= 

1,1 

OSKRD 

=  0 

IBTABL 

= 

1 

NFL  AM 

•  1 

NPTABL 

= 

1 

NEROS 

=  0 

NDYN 

= 

1 

NHTU 

=  1 

NBF  = 

'i  ,  1' 

1 

a  IIRES(  1  :  = 

NRESiat  -  0  Snot  used 

LDBEC 

1 

*  0 

=  PROP 

BED  TNITIALLY  UNCOMPACT,  l=rOMPACTEr 

JHTW 

= 

!  ) 

*1 

LYER 

SIBPES  =  1 

.Y>n  ved  *'c  resi 

NTC 

0 

*  I 

NHie  - 

0.0,0  NC  longer  NEEDED 

NXCU 

= 

0 

NBLOUN 

TS 

0 

NSUSOL 

S 

0 

KMODE 

= 

0 

GUN  aIIqms  modeling  of  carter  hole  i  i" 

MODET 

3 

0 

*  1 

=  TANt 

NECHO 

= 

(1 

INBCX 

=- 

0 

» 1-1  crop 


SXINTEGRA  *M3  INTEGRATION  PARAMETERS 


IDIM 

ITPRT 

15  s4.:. 

0 

•»ao  s 

10 

*20 

ISTEP 

3 

100 

TSTOP 

s 

0 . 02  ' 

lOTST 

3 

0 

ZSTOP 

s 

60.184 

ISTOP 

3 

JSC'O 

TINT 

s 

I’.OOOa 

lUNIN 

« 

0 

SAFE 

3 

1.1  » 

lUNOUr 

CRIT 

R20LV 

j  ,  c 

0 . 05 

*1.1 

*0 .  Ot 


’’ABLIB  =  0.0002 
TABLP  =  .  OOC  1 


*M4  FILE  COUNTERS 


NSTA 

NTEM 

NEPS 

NZPT 

MORE 

NBL 


=  4  moved  to  GUN  SNUMBER  OF  TUBE  DIAMETER  STATIONS 

=  0  NEL  =  O 

-  0.0,0 

=  3  *  MOVED  TO  RES I  DEC! 

=  0,0.0 


*GUl'l 

NAME 

% 

* 

« 

MGAtl 
NSTA 
ZA  = 
SIJND 
TUST 


CWAM 

land 

TWST 

CLEN 


10.0384 
=  1.1311 

=  ^9FP9 . R 

=  7,8^40 


3.4<il6.  25.2362,  51.  .173 


sRESI 


=  'SOMM  SMOOTHBORE’  t 
GRVE  =  1.1311 

G/L  =  l.OOOC 

TRAV  =  53.7P08 

=  5  GLOC  =  -0.0^2®.  -2.0079, 

*  u  *NUM  nf  Stations  in  tube  profile 
.  O  .  "  .  55  ,  2  .  <  '5  ,  60.1  04 

lA  =1.26.  1.26.  1.1811.  1.1811  *  or  use  FA  iradius) 

ss  90990  s  p[r  =  4,.0  3r  use  'LST  (  ri.;ST=360  R  IF  >  1=  c-ns  nght'^T-’T 


tCl2  POPE  FE5'S'"ANCE  -  3^3  =  DICT.  BR  =  PRErSiJRE 
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rSRES  »  I  •LINE'^P  RE5 1  stance  .  a-BA9ED  ‘  ELOCITV 
t  .f  TRAV  or  BR2  it*rt  with  O.C.th#  '-e*t  of  ’•Me  values  are  ■'elstiv* 
«  to  initial  proj.  location,  if  '  0.0  all  is  reiati-e  tc  Pre*?cM 


S  AIR 
NPT5 
TRAV 
PRE3 


■  0.0,  •'.3169,  0.6169,  53 .  ^OOB 

=  O.C'.  1300.0.  100.00.  50.0-M8 


»PR0J  *ni3  PROJECTILE  AND  RIFLING  CHARACTERISTICS 

2BPR  «  8.03  »7,a3  «a.03  *9.836 

PRUT  =  0.723  *  UTPR  or  PRUT  or  WT 

PRIN  »  0.14  *«P0LAR  MOM  OF  INERTIA  --LBM-IM**g 
SRIF*  6.0  movea  '■o  S6UN 


ixahbgas  *m5  amb  jAS  props 

■"■"^■^TEMST  =  320  PST 

IgmST  -  39  GAMST 


-  14.7 


XPROPBED  *M6  PROPELLANT  BED  GENERAL  PROPS 

n*  can  use  either:  TEMPR  .or  temPF.  or,  I'EMPC 
TENPC  -  ai.o 


*PRIM 


■«  SC8  IGNITER  THERMCCHEMtSTPv 

NAME  *  LX  capillary •  S  CHUT  =  .3E046E-O3 

GAMA  »  1.1100  COV  =  33.9743 

TEMP  a  300.00  FQRC  =  23433.3 

*  EIG  3  6303000.  ar  if  FQRC  is  used  it  will  be  'j.sed  internal  Iv 

*  to  calculate  EIG  »  <F0RC*13.D0> . (GAMA-i .DO) 

6MIS  «  36.13 

(NOTE:  Following  is  a  guess  and  mav  not  match  IBHVG 
JJ  *  3  SNUMBER  OF  TIMES 
II  «  3  *NUMBER  OF  POSITIONS  <SEE  M4 i 
[POSIIG  *  0.0,  3.0,  4.00 

ItimEIG  =  0.006.  U.008,  O.OlO 

I  RATEIG  «  0.070,  O.070.  0.070, 

0 . 040 .  0 . 040 ,  0 . 040 , 

0.010,  0,010,  C.OlO  *  •'.a  cpntT.ents  above  this  alloi 


SFROP 


NAME 

= 

•LX  CAPII 

CHUT 

M 

0.0331 

DIAM 

9 

0.5116 

«  WEB 

S 

0.163 

GAMA 

m 

1 .1100 

COV 

9 

33.9743 

TEMP 

s 

300 . 00 

ALPH 
*  NTBL 

=  1 

1.0 

o  comtT.ents 

above  this  allow 

GRAN 

2 

ipc  *2  *' IPF • 

LEN 

= 

3 . 0000 

DP 

* 

0 . 1 873 

RHO 

0.0434 

FCRC 

33455.3 

BETA 

■= 

55.1181 

I^N^e:  'nova  input  needs  follow  for  prop  del 
■•■ISGRLFT  «  0.0  SLEFT  HAND  LOC  OF  BAG  OF  PROP  I 
IZGRRHT  »  3.0  •  is  this  right  ?? 

■RGRl  *  0.09375  SINNER  RAD  OF  PROP  note  this  is  one  dimensional 

■RGRO  *  0.3358  SGUTER  RADIUS  notes  NOT  used  in  cooe 


SLW  =  0.0  SSLOT  WIDTH 

NFIX  =  0  S  0»GRAIN  FREE  TO  MOVE,  1  =  ATTACH  TC 
BONDX  =0.0  *  STRENGTH  OF  BOND 


,  3=  ATT  PROJ 


GAP  =  17400.  SRHEOLOGY  RATE  OF  PROPAGA  OF  STRESS  OF  SETTLED  BED 

GEPQ  ='■>.8  *0.35  *  0.830.38  *0.35  *0.301  *0.58387 

GCAP  =  SO'-'OO.  »RATE  OF  »ROPAGA 

ANU  =  0  *POISSaNS  RATIO  OF  PROP 
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TEMPIG  *  910.  •  IGNITION  TEMPEPATUPE  (3^ 

KP  =  .OS'’"’  *  thermal  conduct  I 'VI  TN 

ALPHAP  »  .0001345  *  THERMAL  DIFPUSIVITV 
EMMIS  -  -to  *  EMMISIVlTv  FACTOR 

♦PROP 

NAME  ■  LX  MOD  TUBE  ' 

CHWT  =  O . 003307 

LEN  ’  7.3441 

gran  «  1  PF  »  =  *  ‘  1  PF  • 

DIAM  =  0.19625  DPERF  =  0.1975 

»  WEB  s  0.008736 

GAMA  *  l.UOO 

COV  =  22.9743  RHO  =  0.0434 

TEMP  *  300. 'O’?  FORC  =  33455.3 

ALPHA  a  0.0  BETA  a  0.1 969 

NTBL  =  1 

t  Mote:  ncva  input  needs  follow  for  prop  ceh 


2GRLFT  X  0.0  SLEFT  HAND  LOC  OF  BAG  OF  PROP  I 

ZGRRHI  X  7.2+413  SRICHT  HAM  BOUND 

R6RI  =  0.1975  SINNER  RAD  OF  PROP  note  this  is  one  dimensional 

RGRQ  =  0.19625  SOUTER  RADIUS  note:  NOT  used  in  code 

SLW  X  0.0  SSLCT  UIDTH 

NFIX  =0  S  OxGRAIN  FREE  TO  MOVE.  1  =  ATTACH  TO  TUBE,  2=  ATT  PRQJ 
BONDX  X  0.0  S  STRENGTH  OF  BOND 

GAP  X  17400.  SRHEOLOGY  RATE  OF  PPOPAGA  OF  STRESS  OF  SETTLED  BED 

GEPQ  xO.a  SO. 35  .  S  0.920.28  SO. 25  SO. 301  SO. 58387 

GCAP  X  50000.  SPATE  OF  PROPAGA 

ANU  X  0  SPOISSONS  RATIO  OF  PROP 

B1  X  0.0  s  BURN  RATE  CONSTANT 

TEMPIG  *  910.  S  IGNITION  TEMPERATURE 

KP  X  .0277  S  THERMAL  CONDUCTIVITY 

ALPHAP  »  .0001345  S  THERMAL  DIFFUSIVITV 
EMMIS  X  .6  S  EMMISIVITY  FACTOR 


SPROP 

NAME 

X  ' wcass ' 

‘oRAI'l 

X  BALL  S18 

S  IGNC 

X  1 

THRC 

*  .5000CE-O3 

NTBL 

X  1  S-1 

DIAM 

*  0.0630 

S  DL/S  X  0.008570,  0.0315,  0.0851  '  NOVA  does  not  like  's 

DLS  X  0.008570,  0,0315,  0.0851 

S  DEPL  =  0.0005399,  0.0019845,  .0053613 

RHOL  X  0.0539.  0.0539,  0.0548,  0.0589 

GAML  X  1.2850,  1.2850,  1.2790,  1.2340 

FRCL  X  254259.?,  254259.9,  278213.9,  373962.8 

COVL  *  30.8077.  30.8077,  30.9184,  26.4897 

TMPL  X  1896.=.  1896.9,  2109.4,  3471.2 

COEL  X  .62236E-03,  .62236E-03,  .75113E-3,  0.001788 
eXPL  X  0.8053,  .8053,  .8053,  .8053 

CHWT  a  0 . 3336 

S  Note:  nova  input  reeds  fellow  for  prop  dek 

ZGRLFT  X  >,0  ILEFT  HAND  LOC  OF  BAG  OF  PROP  I 
2CRRHT  X  7.5  Sguess  SRIGHT  HAN  BOUND 

RGRl  X  0.25  ScuessSINNER  RAD  OF  PROP  note  this  is  one  dimensional 

RGRO  X  0.63  sCLTER  RADIUS  note:  NOT  used  in  code 

SLW  X  0.0  »SLCT  WIDTH 

NFIX  X  0  S  0=GRAIM  FREE  TO  MOVE.  I  *  ATTACH  TO  TUBE,  2*  ATT  POd 
BONDX  =  0,0  S  3TPENGTH  OF  BCND 

CAP  X  17400.  SRHEOLOGV  RATE  OF  PROPAGA  OF  STRESS  OF  SETTLED  BED 
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acv  .  — v; 


5CAF  *  30000  JRATE  OF  »ROPAGA 

»NU  *  0  •pnrssoMC  rat  re  of  prof 


@ 


=  0.:-  *  3UPN  FATE  CONSTANT 

rEMPIG  =  810.  «  rGN!T10|.i  TEMPERATURE 

F  »  .  02“'"  t  THERiiAL  COMDL'C’ I V I T  / 

^LFHAP  »  .O'Olri-S  »  rH£PMA(_  D  IFFiJ^  T.  I  T  < 

tMMI3  ^  .3  •  EMMISIV.T/  FAi-TGR 


fCOMM 

i  FOUR  LA/ER  PROPELLANT  iODEL  iNONINAL  1  MM  T/IAM* 
i  warning  this  dek  is  .nec-ic  because  was  ramerteo  cut  durirn  cen/ers 
JAME  *  • UC685 ‘  gram  =  BAL  L 

IGNC  =  1  TWRC  =  O.OOOS 

DIAM  «  0.00100838 

NTBL  a  -I 

DL/S  a  0.01 3fc,  O.OSOO,  0.1330 

RHQL  a  1498.6,  1498.6,  1513.8.  1631.3 

GAML  a  i.BBS,  1.885,  1.879,  1.834 

FRCL  a  760000,  760000,  831600,  1117800 

CQVL  a  0.001113,  0.001113,  0.001  ll'T.  0.000957 

TdPL  a  1896.9,  1996.9,  8109.4.  3471.2 


CFIL 

a  0 . 00087 

EXIL 

a  0.8053 

CFSL 

a  0.00087 

EX8L 

a  0.8053 

CF3L 

a  0.00105 

EX3L 

a  0.8053 

CF4L 

a  0.00850 

EX4L 

a  0.8053 

CHWT 

a  0.1513 

♦ETC 

♦  warning  'his  is  metric  et  dei.  ****»*»♦***♦*♦******♦•**•*••*** 

PFACal 
NPWP  a  9 

TPWR  a  O.  3E-5,  6E-5.  1.8E-4,  lE-3,  1.85E-3,  1.38E-3,  l.^E-a.  8E-5 
PWR  a  O.  0,  4.4E6,  8.4E6,  4.6E6,  1.83E6.  0.43E6.  0.09E6.  ' 

♦END 
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CHALLENGES  FOR  ETC  MODELING 
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ET('  Gun  Modeling  at  CRT: 
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Unitlixedness  Modeling  Using  PDF  Approach  (continued) 


Turbulent  Unmixedness  Modeling: 
Choice  of  States  for  a  Two  Species  System 
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Solution  Procedure  for  Probabilities 
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Solve  for  Lagrange  multipliers  using  Newton's  method,  e.g. 


Temperature  Spottiness  Modeling 


i.e.,  maximize  the  entropy  of  mixing.  Again,  using  Lagrange  multipliers,  solve  for 


mperature  Spottiness  Modeling  (continued) 
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ETC  Gun  Modeling  Approach 


Entrainment  of  droplets  from  working  fluid 


ETC  Gun  Modeling  Approach 

(continued) 
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ETC  Gun  Modeling  Approach 

Schematic 
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ETC  Gun  Calculation  Description 
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MAZe  Code  Simulalfcm  of  30mm  ETC  Gun  Test 

Selected  Densities  at  0.50  msec 

Turbulent  Unmixedness  Case 


Distance  (mm) 


^ode  Simulation  of  30mm  ETC  Gun  Test 
Droplet  Density  Time  =  0.50  msec 


Distance  (mm) 


MAZe  Code  Simulation  of  30mm  ETC  Gun  Test 

Droplet  Extinction  Time  Time  =  0.50  msec 
Turbulent  Unmlxedness  Case 


r 
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Distance  (mm) 


M  AZe  Code  Simulation  of  30mm  ETC  Gun  Test 
Turbulent  Damkohler  Number  (W^etno)  msec 


Distance  (mm) 


Figure  17.  Comparison  of  luriwlent  unmixedness  and  perfect  mixing  cases  at 
1.1  msec:  total  fluid  density  and  product  density  contours. 


TURBULENT  UmilXEDNCSS  CASE 


Figure  6.  Projectile  base  pressure  and  acceleration  time  history. 


MAZe  Code  Simulation  of  30mm  ETC  Gun  Test 
Comparison  with  Experimental  Overpressure  Time-Histories 
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Time  (msec) 


MAZe  Code  Simulation  of  30mm  ETC  Gun  Test 
Comparison  with  Experimental  Overpressure  Time>Histories 
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1 

Princeton  Combustion  Rsch  Lab 
Princeton  Corporate  Plaza 

ATTN:  Dr.  Neale  A  Messina 

11  Deeipark  Dr..  Bldg.  IV,  Suite  1 19 
FAmmouth  Junction,  NJ  0^52 

1 

Defenae  Nuclear  Agency 

ATTN:  COL  Downie 

6801  Telegirqjh  Rd. 

Alexandria.  VA  22310 

1 

Science  Applications  IntT  Cotp. 
ATTN:  Dr.  Gary  Philip 

10210  Campus  Point  Dr. 

San  Diego,  CA  92121 

1 

Director 

Sandia  National  Laboratories 

Computational  Fluid  Dynamics 

Div.  1511 

ATTN:  Dr.  Gr^oty  F.  Homicz 

1 

Science  ^jpUcations  IntT  Cotp. 
ATTN:  Dr.  C.-C.  Hsiao 

10210  Canqws  Point  Dr. 

San  Diego,  CA  92121 

Albuquerque,  NM  87185 

1 

Science  Applications  IntT  Coip. 
ATTN:  Dr.  RedSu 

1 

Director 

Sandia  Natkinal  Laboratories 

Advanced  lYojects  V 

10210  Campus  Point  Dr. 

San  Diego.  CA  92121 

ATIN:  Dr.  Steve  Kempka 

Albuquerque.  NM  87185-5800 

1 

Science  ^jplications  IntT  Cotp. 
ATTN:  Dr.  JadBatteh 

1519  Johnson  Ferry  Rd.,  Suite  300 

1 

Defense  Nuclear  Agency 

ATTN:  Dr.  R.  Rohr 

Marietta.  GA  30058 

6801  Telegraph  Rd. 

Alexandria,  VA  22310 

1 

Science  ^tplicadons  IntT  Corp. 
ATTN:  Dr.  Sanford  Dash 

501  Office  Center  Dr. 

1 

Hie  PBnns]rivania  State  Uraversity 

ATTN:  Fnrfessor  F.  B.  Cheung 

Prop  Engr  Rsdi  Ctr/Mech  Engr  Dqx 

Suite  420 

Fort  Washington.  PA  19034-3211 

140  Research  Building  ”E” 

Universiqr  Rsrfc.  PA  16802 

1 

Science  Applkatioas  IntT  Coq>. 
ATTN:  Mr.  Neerqi  Sinha 

501  Office  Center  Dr..  Suite  420 

1 

Universiqr  of  Tennessee  Space  Institute 

Center  for  Laser  Applications,  MS-14 

Fort  Washington,  PA  19034-3211 

ATTN:  Dr.  Dennis  Keefer 

TuOidioma.  TN  37388-8897 

1 

Science  ^jplications  IntT  Cotp. 
ATTN:  Dr.  Ashwin  Ifasangadi 

501  Office  Center  Dr. 

1 

The  Pennsjrivania  State  Umversity 

Prop  Eogr  Rsch  Ctr/Mech  Engr  Dqtt 

ATTN:  Professor  Kenneth  K.  Kuo 

SnilB420 

Fort  Wadungton,  PA  19034-3211 

140  Research  Buildii«  "E" 

University  Park,  PA  16802 

1 

FMC  Corporation 

Naval  Systems  Division  -  M170 
ATTN:  Mr.  RMrick  Janke 

4800  E  River  Rd. 

Minneapolis.  MN  55421-1498 

29S 


Naof 


No.  of 


COwes  Omwiatioii 

1  FNK:  Coqxmion 

Nival  SyAems  Diviaioa  -  M170 
ATIN:  Mr.  Jahn  Dyvik 
4800  E.  River  Rd. 

Miniieapolis.  MN  S542M498 

1  Techncdogy  Gateways  Ricaqxxated 
ATTN:  Mr.  Jete  Brown 
600  East  John  Sims  Parkway,  Suite  118 
NioeviUe.FL  32S78 

1  Veritay  Technology  Inc. 

ATTN:  hfr.  Edward  B.  Fisher 
484S  Millerqiort  Highway 
P.O.  Box  30S 

East  Amherst,  NY  14031-030S 

1  GT-Devices,  Inc. 

ATTN:  Dr.  Niels  K.  Winsor 
S70SA  General  Washington  Dr. 
Alexandria,  VA  22312 

1  GT-Devices,  Inc. 

ATTN:  Mr.  Dave  Horiburt 
S057A  General  Washingtoa  Dr. 
Alexandria,  VA  22312 

1  GT-Devices.  Bic. 

ATTN:  JonEarahan 
S70SA  General  Washingion  Dr. 
Alexandria.  VA  22312 

1  GT-Devices,  Inc. 

ATTN:  Dr.  J.  Robert  Greig 
S70SA  General  Washington  Dr. 
Alexandria.  VA  22312 

1  M  Rteedman  ft  Associates 
ATTN:  Dr.  Eli  Freedman 
2411  Diana  Rd. 

BaltimQte.MD  21209-1S2S 

1  Paul  Gough  Associates.  Inc. 

ATTN:  Dr.  Paul  S.  Gough 
1048  South  St 
POrtamouth,  NH  03801 


Copies  Organization 

1  California  Research  ft  Technology  Division 
The  TITAN  Corporation 
ATTN:  Dr.  Phil^  A.  Hookham 
20943  Devonshire  St 
Chatsworth,  CA  91311-2376 

1  Olin  Ordnance 

ATTN:  Mr.  Karl  R.  Linde 
200  E.  High  St 
Red  Lion.  PA  173S6 

1  Olin  Ordnance 

ATTN:  Mr.  Hugh  A.  McElroy 
10101  9th  St  North 
St  Petersburg,  FL  33716 

1  State  University  of  New  Yoric  at  Buffalo 

Dq;>artmem  of  Electrical  Engineering 
ATTN:  Professor,  Dr.  W.  J.  Saijeant 
312  Bonner  ECE-SUNY/AB 
Buffalo.  NY  14260 

1  General  Dynamics 

Land  Systems  Division 
ATTN:  Dr.MelWidner 
P.O.  Box  2074 
Warren.  MI  48090-2074 
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Dir.USARL 

ATTO:  AMSRL-WT-P.  I.  M«y 
AMSRL-WT-PA. 

T.  Coffee 
G.  Wien  (IS  q>s) 

P.  Tran 

S.  Richardson 
W.Oberle 

K.  White 
A.  Juhasz 
AMSRL-WT-PB. 

E.  Schmidt 

C.  Nietubicz 

AMSRL-WT-PC.  R.  Fifer 
AMRSL-WT-PD.  B.  Bums 
AMSRL-WT-PE. 

M.  Nusca 

F.  RoU^ 

P.  Conroy 

T.  Minor 

D.  Kooker 
O.Keiler 
A.Haat 

AMSRL-WT,  D.  EockAaU 
AMSRL-WT-T,  W.  Monisoo 
AMSRL-WT-WD, 

J.PoweU 

A.Niiler 


B 


iNIEhmONALLY  LEFT  BLANK. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undeitakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
commentsAuiswers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ABL-SR-IO _ of  Report  May  1994 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ ^ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Irrdicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City.  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correctirai,  please  provide  ttie  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  uqre  dosed,  and  mail) 
(DO  NOT  STAPLE) 


